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The Jarzynski relation is a recently discovered result relating the average exponential of the work 
done under nonequilibrium conditions to an equilibrium free energy difference. We illustrate this 
remarkable relation by considering the expansion and compression of a classical ideal gas particle 
inside a cavity with a moving piston in one dimension. Using a simple relationship between the 
distributions of the nonequilibrium work done in the two (forward and reverse) processes, we verify 
the Jarzynski relation and discuss the practical realization of the Jarzynski equality for the model. 



I. INTRODUCTION 

Most relations in thermodynamics are valid under 
equilibrium or near-equilibrium conditions. For exam- 
ple, consider a system in contact with a heat bath at 
temperature T. If under some (infinitesimal) process an 
energy Q is exchanged with the bath, the first law of 
thermodynamics tells us that 



Jarzynski announced to the scientific community the 
following equality that now bears his name 



Q^AU + W 



(1) 



where AU is the change in the internal energy of the 
system and W is the work done by the system. The 
statement is of course an expression of the conser- 
vation of energy, which is valid not just in equilibrium 
or whether the process is reversible or irreversible. If we 
assume that the process takes place very slowly (qua- 
sistatic, reversible), the second law of thermodynamics 
further relates the entropy change AS to the energy Q 
via the relation Q — TAS. With this replacement, the 
relationship may be recast as 

-W = AU - TAS = AF (2) 

relating the work done on the system — W to the change 
in the Helmholtz free energy AF. 

If the process is not quasistatic, we have from the sec- 
ond law of thermodynamics, (Q) < TAS, and the equal- 
ity ^ converts into an inequality 



(W) > AF 



(3) 



where the < • • • > denotes an average over an ensemble 
of measurements of W (or Q), where each measurement 
is taken as the system undergoes the same variation in 
external parameters starting from an equilibrium state. 
(For the system of a gas under a piston and the process 
consisting of movement of the piston, this means that 
for each measurement the piston moves along the same 
trajectory, with the gas in thermal equilibrium at the 
start of every trajectory. The external parameter in this 
case is the volume of the enclosure, which is the same at 
the start of every trajectory and the same at the end of 
every trajectory.) 

Many macroscopic physical quantities measured in ex- 
periments where the system is driven away from equilib- 
rium satisfy inequalities such as (jJJ). However, in 1997, 



,W/kBT 



(4) 



where the angle brackets have the same meaning as in 
the previous paragraph. This remarkable formula tells 
us that we can extract equilibrium information (AF) 
from the ensemble of nonequilibrium (finite-time) mea- 
surements Pj. By application of the mathematical iden- 
tity (e^) > e^^^ the Jarzynski relation Q yields the in- 
equality in Q. 

For our purposes, the Jarzynski relation is valid when 
at the start of every measurement the system is at ther- 
mal equilibrium with a heat bath, and that during the 
finite-time switching of the external parameters, no en- 
ergy is exchanged with the heat bath (i.e. the process 
is adiabatic). The free energies Fi and Fq that comprise 
the difference AF — Fi — Fq refer to states of the system 
in equilibrium at the external parameters corresponding 
to the start and end of the switching process. More de- 
tailed and general conditions for the validity of I^J can 
be found in 0,111111 [13. 

The utility of the Jarzynski relation lies in its potential 
to give a good estimate of the free energy change, AF, 
given a finite number of measurements of the work within 
a finite amount of time. (e.g. compared with the use 
of 0) In fact, the Jarzynski relation has been used to 
measure free energy differences in experiments involving 
single molecule systems such as the pulling of RNA or 
protein [lHllla. 

Nevertheless, it should be pointed out that there ex- 
ists considerable recent debate in the literature concern- 
ing when precisely the Jarzynski relation is valid and 
whether it is applicable in practice, even for an ideal gas 
E in 111 III 111 El- In reference [l^, Lua and Gros- 
berg discussed the practical applicability of the Jarzynski 
relation in the context of an ideal gas particle expanding 
under a piston and moving in one dimension. (See figure 
n The Jarzynski equality for N independent particles is 
given by the Jarzynski equality for 1 particle raised to the 
power N.) In that work, the Jarzynski relation was ver- 
ified with an explicit calculation. The work distribution 
function was also derived and given explicitly. Finally, 
the limit of a fast moving piston was examined, leading 
to the conclusion that it might be difficult to realize the 
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Jarzynski equality in practice for the case of an expan- 
sion with the fast moving piston, as was also pointed out 
in 0. 

Our goal in this article is to illustrate the Jarzynski 
relation and gain some insight into how it work s by con- 
sidering the same ideal gas model discussed in [lOj . Our 
results and conclusions for the case of an expanding cav- 
ity are the same as in that work. However, in the present 
work we also examine the 'reverse' case in which the pis- 
ton compresses the cavity. Furthermore, we verify the 
Jarzynski relation in an indirect but simpler fashion. 

The plan of the rest of the article is as follows. In sec- 
tion ^ we give explicit expressions for the probability 
distribution functions of the work done by the gas parti- 
cle (or piston) for the case when the piston expands the 
cavity (arbitrarily called the 'forward' process, with dis- 
tribution denoted by PpiW), where the argument W is 
the magnitude of the work done by the gas) and for the 
case when the piston compresses the cavity (called the 
'reverse' process, with distribution denoted by Pr{W), 
where W is the work done by the piston). These proba- 
bility distributions appear when calculating the average 
defined above, e.g. {e^/''^^)^, = PFiW)e^^''"^dW. 
We assume that the piston moves with the same constant 
speed during the expansion and compression. In an at- 
tempt to make the main presentation short and simple, 
we shall also relegate cumbersome derivations to the Ap- 
pendix. In section ITTTl we show that the forward and re- 
verse distributions, Pf{W) and Pfi{W), satisfy a simple 
relation, which we then use to show that the Jarzynski 
relation is valid for our model. In section Hvl we discuss 
the practical evaluation of the free energy difference from 
the exponential average in the Jarzynski relation for the 
forward and reverse processes. 



II. PROBABILITY DISTRIBUTION OF THE 
WORK 

Referring to figure ^ the initial state of the particle 
(x, v) consists of a position x distributed with uniform 
probability within the initial length of the cavity, and a 
velocity v drawn from a Maxwell-Boltzmann distribution 
(i.e. a Gaussian distribution, e"" /^). To simplify the 
writing and without loss of generality, we use units such 
that ksT = 1, the mass of the particle to = 1, and that 
the piston displaces by an amount AL during time r = 1. 

For the forward process in which the cavity expands 
from 'volume' L to a final volume of L + AL, the pre- 
scription for evaluating the distribution is 



Pf{W) = 



1 



2nL 



dx 



dve-''^/^5{W -Wr{x,v)) 



(5) 

where Wr{x, v) is the work done by the particle given the 
initial coordinate {x,v). The 6 function picks out the 
initial states which lead to an amount of work W. The 
calculation is presented in detail in Appendix ^ The 



final result is 



Pp{W) = S{W)Po + ^-^^ f{W) . (6) 



2TTnAL 

Here, Pq is the probability to obtain vanishing work be- 
cause the particle is unable to chase the piston or hit it 
even once. 



1 



dx 



(L+AL) 



rf„e-(---)'/2 . (7) 



^/QjtL Jo J-(l+al) 

The integer n (which is the number of bounces by the 
particle against the piston) is obtained in Appendix ^ 
and is given by the formula 




2W 



AL(2L + AL) 




(8) 



where [. . .] means integer part of . . .. For example, simple 
algebra indicates that as long as < 4AL(2L -I- AL), 
we have just one collision, n = 1. For the values of 
work W in the next interval, 4AL(2L -I- AL) < W < 
12AL(2L -I- AL), we have n = 2, etc. 

Finally, the function f{W), which we call the overlap 
factor, varies between and 1 and is illustrated in figure 
121 An explicit expression is given by ljA10|l in Appendix 

El 

Similarly, for the reverse process in which the cavity 
contracts from 'volume' L + AL to a final volume of L, 
the prescription for evaluating the distribution is 



1 



L+AL 



dx X 



27r(L + AL) Jo 

dve-'^'/^SiW -w'^{x,v)) , (9) 



Pb.{W) - 



where 'w'^{x,v) is the work done by the piston given the 
initial coordinate of the particle {x,v). The final result 
is 



Pr{W) = s{w)p;, 



Here, Pq is given by 



■l(nAL- 



Lf{W) 



27mAL (L + AL) 



(10) 



P' 

^0 



1 



27r(L + AL) 



L+AL 



dx 



(11) 

The integer n and the function f{W) are the same as in 
the forward case. 

Thus, the probability distributions ((HJ and H10|l consist 
of a (5 function peak at = and a tail or hump at pos- 
itive W (see figureOfor the case of a fast moving piston). 
The forward distribution is examined in Appendix [21 in 
the limit of a fast moving piston and in Appendix ^ in 
the limit of a slow moving piston. 
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III. VERIFICATION OF THE JARZYNSKI 
RELATION 

We now show that our model satisfies the Jarzynski 
relation Taking the ratio of expressions © and (|10|l . 
we obtain 



Pf{W) 
Pr{W) 



(12) 



where AF 



ln( 



is the difference between the 



free energies at volumes L + AL and L, in equilibrium at 
the same temperature. (Relation (jl2|l is reminiscent of 
another remarkable result called the fluctuation theorem 
or Crooks relation, which relates the distribution of en- 
tropy productions of a driven system that is initially in 
equilibrium to the entropyproduction of the same system 
driven in reverse H i, 11 Q, m . The simple proof below 
was also inspired by the first two of these references.). 

The exponential average of the work is readily evalu- 
ated from (|12|l and the normalization of the probability 
distributions. For the forward process, 

Jo 

PR(W)dW 



10 

-AF 



-AF 



(13) 



Similarly, for the reverse process 
Jo 

PF{W)dW 



/•OO 

10 

„AF 



„AF 



(14) 



The Jarzynski relation is also verified explicitly in ap- 
pendices [21 and ^ in the opposite limits of a fast moving 
piston and a slow moving piston. 



IV. REALIZING THE JARZYNSKI EQUALITY 
IN PRACTICE 



In this section we examine how well the average of 
the left-hand-side of the Jarzynski relation (@J) can repro- 
duce the equality with the right-hand-side, based on a 
finite-number of measurements of the work. This issue 
is of immense practical interest in experiments involv- 
ing the manipulation of nanoscale objects and biological 
molecules. 

From the definition of the exponential average, 
(e^)^ = PFiW)e^dW, one can already see the po- 
tential problem. In an expansion, when the piston moves 
very fast, only the very fast (and very rare) particles can 
hit the piston to produce non-zero work. But these fast 



particles can give a significant contribution to the aver- 
age because of the factor e^. Many measurements or 
repetitions of the expansion are necessary to sample an 
adequate number of particles that move fast enough in 
order to assess this contribution accurately. 

To be more precise, and to see a clear difference be- 
tween the forward and reverse processes, we consider the 
case when the piston moves very fast, such that a typical 
particle moves very little during the expansion or com- 
pression of the cavity. This is equivalent to the conditions 
L ^ CT and AL ^ cr, where c = y/kBT/m is roughly 
the speed of sound in the gas. In our simplified units, 
this translates to L, AL ^ 1. 

First, let us rewrite H12|) as 



-AF 



(15) 



we see that the contribution of the expanding or forward 



work measurements to the average, 



-AF 



is given by the 



reverse work distribution, Pii{W). This indicates that 
measurements of Pf{W) can give an accurate result for 
the average if Pf{W) is significant at values of W where 
Pr{W) is significant; that is, if there is a reasonable over- 
lap between the forward and reverse work distributions 

m 

Now let us determine the probability of obtaining zero 
work for the forward and reverse processes. For the for- 
ward process, the probability -Pf(O) is clearly close to 1. 
From the relation H12|l . the probability of obtaining zero 
work for the reverse process is given by. 



Pb{0) 



L 



AL 



(16) 



Clearly, in the limit of a fast moving piston, this proba- 
bility is given by the fraction of the cavity volume that 
is not swept or visited by the piston. 

As to the probability of obtaining non-zero work, in 
the forward process this clearly vanishes (the area under 
the lower curve in figure ^ in the limit of a fast moving 
piston. In contrast, the probability of obtaining non-zero 
work in the reverse process is finite (the area under the 
taller curve in figure This probability is given by the 
fraction of the cavity volume that is swept by the piston, 
AL/{L + AL). 

Now we can address the issue presented at the start 
of this section. In the forward process, rare events (non- 
zero work values) contribute a finite fraction to the aver- 
age, given by the area under the 'reverse' curve of figure 
|31 This means that many measurements are necessary 
(given roughly by the inverse probability of obtaining 
non-zero work values ljC3|l ) in order to faithfully produce 
the equality in the Jarzynski relation I^J. In contrast, for 
the reverse process, the dominant contribution is given 
by the zero work values, which have finite probability. 
Therefore a smaller number of measurements in the re- 
verse process can reproduce the equality in (0J quite well, 
as was indeed revealed in computer simulations. 
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V. CONCLUSION 

We obtained explicit expressions for the distribution 
of the work done in the forward and reverse processes, 
Pp{W) and P[^{W), for a simple ideal gas model. By 
noting a simple relationship between the two distribu- 
tions (equation H12|l ) we illustrated the validity of the 
Jarzynski relation (equation for the model. We also 
conclude that in the forward or expanding case, the eval- 
uation of the average exponential of the nonequilibrium 
work from a finite number of trials or measurements can 
give a poor result for the free energy difference when the 
piston moves sufficiently fast (indicating that repeated, 
rapid nonequilibrium measurements may not be auto- 
matically advantageous [13 )• This is in sharp contrast 
to the average that can be obtained from the comple- 
mentary or reverse process. In general an optimal result 
may be obtained by a suitable combination of the mea- 
surements in the forward and reverse processes |17| . 
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APPENDIX A: CALCULATION OF Pf{W) 

In this section, we first derive an expression for Pp{W), 
based on the work in llClj |. Again, we take reduced units 
such that ksT = 1, the particle mass to = 1, and the 
piston moves with speed AL/r during a time interval 

T=l. 

Let us first assume a positive initial velocity (refer to 
figure^, in which the particle can strike the piston first 
before hitting the left end of the cavity. The time taken 
for the first collision with the piston is ti — J"^^r ■ After 



the collision, the velocity of the particle relative to the 
piston gets reversed and the speed of the particle gets 
diminished to u — 2AL (assuming v > 2AL). The time 
taken for the second collision with the piston is given by 



t2 — t,'^3AL • general, for the n*'' collision 



c = 



(2n- l)L-x 



(AI) 



(2n- 1)AL 

Similarly, for a particle with a negative initial velocity, 
(2n- l)L + x 



v-{2n- 1)AL 



(A2) 



These relations can be inverted to give conditions that 
should be satisfied by the speed of the particle in order 
to result in exactly n collisions with the piston within a 
time interval r = 1. For positive initial velocities, 

{2n-l){L + AL)-x < \v\ < {2n + l){L + AL)-x (A3) 

For negative initial velocities, 

{2n-l){L + AL) + x < \v\ < {2n+l){L + AL)+x (A4) 

The work done by the piston on the particle after one 
collision is the change in momentum of the particle times 
the velocity of the piston, 

-wi = {-{v - 2AL) - v)AL = -2{v - AL)AL (A5) 

In general, the work done after n collisions is 

-Wn = -2vALn + 2AL^n^ (A6) 

Note that the work done can also be calculated from the 
change in kinetic energy after n collisions, 

-Wn = ]^{v-2nALf-v^ /2 = -2nALv+2n^ AL^ (A7) 

The work done by the particle on the piston is positive 
for an expanding volume. 

The inequalities (|A3|) and IIA4|l lead to the following 
partition of the integral in 



J 



Pf{W) = 



^ da; V / dwe-^'/2 ^ s {W - {2vALn - 2AL^n^)) 

^ da; V / dwe-"'/2 x 5 {W - {2vALn ~ 2AL^n^)) 



dx 



V2ttL Jo J -(l+al)-x 



){L+AL)+x 
(L+AL)-x 



dve-'' /^6{W - 0) 



Call the first term Ii and the second term I2 (the third term is 'trivial'). Performing a change of variable to re- 
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move X from the limits, integrating over x and taking 
advantage of the delta function results in 

h ^ ^ X 



1 



X — X {overlap between 
2/7 



W 



,nAL 



W 



2nAi' ^ ' 2nAL 
and [{2n - 1)(L + AL), {2n + 1)(L + AL)]} 



and similarly for I2, 

00 ^_i(„AL+, 



27rnAi 



E 

n=l 

X ^ X {overlap between 



nAL 



L,nAL' 



W 



2nAL 2nAL 
and [(2n - 1)(L + AL), {2n + 1)(L + AL)]} 

(For two intervals [a, 6] and [c, d] where a < b and c < d, 



the expression {overlap between[a, 6]and[c, d]} equals 
when b < c, equals b — c when a < c < 6 < d, etc.) 
Ji and /2 can be combined as follows, 



h+h = 



- g-i(„AL+^)^ 

L — = ^ 



27rnAi 



X — X {overlap between 



lAL 



W 



- L,nAL 



W 



L 



2nAL ' 2nAL 
and [{2n - 1){L + AL), {2n + 1){L + AL)]} 

00^ -i(„AL+^)^ 

^ V2^nAL ^ ^ 



where the overlap factor f{n,W) satisfies < / < 1, 
since the range of the smaller interval is at most 2L. f 
is also zero for negative W , or positive work values —W 
done by the piston. The conditions that must be satisfied 
by W in order for the overlap associated with integer n 
to occur are 



2nAL {2{n - 1)L + (n - 1)AL) <W < 2nAL (2(n + 1)L + {n + 1) AL) 



(A8) 



Notice that the left boundary of the interval is a function 
of n{n— 1), while the right interval is a function of n{n + 
1). Therefore the right boundary can be transformed 
into the left-boundary by making the replacement n 
n—1. This implies that the intervals HA8|I are contiguous 
and nonoverlapping and that at most one term in the 
summation in Pf{W) survives. One can solve for the 
integer n by taking the integer part (or floor function) of 



a solution to a quadratic equation. 



W > 2nAL {2{n + 1)L + {n + 1)AL) 



(A9) 



which results in formula ijSJ in the main text. 

Dropping the argument n, an explicit expression for 
fiW) is 



J 



fiW) = { 1 



1)( 



AL 
2L 



W 



4nALL 



(" + i)(lr + i)- 



4nALL 



when 

when {n - 1)(AL + 2L) + 2L < 

when {n-l){AL + 2L) + 2L + 2AL < <{ 



{n - 1){AL + 2L) < <{n- 1){AL + 2L) + 2L 

< (n - 1)(AL + 2L) + 2i + 2AL 
^ ' l)(Ai + 2i) 



r 



(AlO) 



APPENDIX B: OBTAINING Pr{W) FROM Pf{W) L + AL. The result is 



To obtain Pi^{W) from Pf{W), make the following 
series of replacements W —W, AL — AL, L 



6 



g-i(-„AL+^) = 

Pr{W) = Pl,5{W) + ^ X 



1 



2TTnAL 



X < overlap between 



-nAL 



2{L + AL) 
W 



2nAL 



- {L + AL),-nAL- 



W 



2nAL 



L + AL 



and[(2n- l)(i), (2n + 1)(L)] 



where 



P' — 
^0 — 



and 



1 nL+AL pL 

/ dx / dve-^^-^^"/' (Bl) 



1 + A 1 



2W 



AL{2L + AL) 




(B2) 



Notice that n on the forward and reverse cases are iden- 
tical. 



The double integral in the Pq term in the forward pro- 
cess is identical to the double integral in the Pq term due 
to reflection symmetry of the integrand, which means 



Po 
P' 



L + AL 
L 



(B3) 



Furthermore the overlaps of the intervals in the for- 
ward and reverse processes are identical. To see this, 
shift each interval in the overlap expression for the re- 
verse case by {2n— 1)AL (this does not change the value 
of the overlap) 



overlap between 



overlap between 



lAL 



~nAL 

W 
2nAL 



W W 
— — -(L + AL), ^nAL + — — 
2nAL ^ 2nAL 



L + AL 



and[(2n- l)(L),(2n + l)(L)] 



- L - 2AL,nAP^ 



W 



2nAL 



+ L 



and [{2n - 1)(P + AL), {2n + 1)(L -f AL) - 2AL] 



(B4) 



Comparing this expression for the overlap in the reverse 
case with that of the forward case, one sees that they 
are identical apart from the shift of — 2Aiy in the left 
boundary of the first interval and a shift of — 2AL on 
the right boundary of the second interval. Because the 
difference of the lengths of the two intervals is 2AL, the 
overlaps are identical. 



APPENDIX C: 



LIMIT OF A FAST MOVING 
PISTON 



In this section, we obtain the form of the forward dis- 
tribution © when the piston takes a very short time r 
to undergo the displacement AL, such that a typical par- 
ticle moves very little during the expansion. That is, we 
assume that AL ^ cr and L ^ cr, where c — k bT Ijn. 
The formulas in this section are similar to that in !Tc|. 

In this limit, the distribution is dominated by a single 
bounce, n = 1, and values W < 4(AL)L, giving 



P{W) ~ d{W)Po 



W 



(CI) 



y/2^AL 4(AL)L 
This form of the distribution suffices to produce the 



Jarzynski equality: 



POO 

/ P{W)e^dW 
Jo 



= Pn 



1 



27rAL Jo 



e sKl2^ ' X 



W 

x——dW 
AALL 



dx 1 
L V27r J- 



V2^AL J- 
X — ^ — dW 



e SAl2 ^ ' X 



= 1 



4(AL)L 
AL 
17 



-AF 



Where AP = — In (1 -f- AL/L). Here, we made approx- 
imations in both the Pq term, by extending the integral 
limits to (— cx3,cx3), and in the tail term, by setting the 
integral lower limit to — oo. 

Let us also calculate the probability of obtaining non- 
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zero work values. Using expression (|C1|I . we have 



Pi 



w>o 



0+ 



P{W)dW 



-AL'^/2 



4\/2^L(AL)2 



W X 



(C2) 



For large AL, we neglect the term W'^/{8AL^) in the 
exponent upon integration, yielding 



Pi 



1 



w>o 



-ALV2 



27ri(Ai)2 



(C3) 



This quantity represents the area under the smaller curve 
in figure 13 



APPENDIX D: LIMIT OF A SLOW MOVING 
PISTON 

In this section, we obtain the form of the forward dis- 
tribution © when the piston takes a very long time r to 
undergo the displacement AL, such that a typical parti- 
cle bounces off the piston many times. This means that 
AL <C CT and L <C cr, where c = ^Jk^frjm. 

Defining AL/L ~ a, we have Pq — and n is very 

large with nAL ~ ^ / Tr-^-v-rxa. Thus, 



2a(2+a) 



P{W) ~ ^ ^ Xtt exp 



aVW 



(1 



a{2 + a) 



w]f{w) . 



(Dl) 

The factor f{W) behaves like a rapidly oscillating func- 
tion in this limit (albeit the period increases as n in- 
creases). We shall replace f{W) by its average value, 
given by the ratio of the area of a trapezoid in figure El 
to the length of the base of the trapezoid, when we 
evaluate the integrals below. 

One can derive the W dependence of this distribu- 
tion using adiabatic invariants, except without the factor 
f{W), as follows. 

In :2j , Jarzynski considers as an example a single par- 
ticle bouncing around inside a three-dimensional cavity 
with hard walls, where the shape of the cavity is a func- 
tion of a control parameter A. In it, he mentions the 

quantity H^^^V\ as an adiabatic invariant (i.e. a constant 
when the control parameter is varied very slowly), where 



H\ is the particle energy and V\ is the volume of the cav- 
ity. These give the work performed on the particle with 
initial energy Eo: ^ Hi - Eo ^ [(Vb/Vi)^/^ - l]Eo 
, and the work distribution when A is switched infinitely 
slowly from to 1: 



lim p{W,ts) 

is — *oo 



1/2 



exp 



(5W 



e(W) (D2) 



where r = {Vo/Vif/^ - 1, /3 = 1/fcsP, and Vb > Vi (the 
cavity is being compressed). 

We would like to employ the same idea here for the 
one-dimensional particle. 

1/2 

In one dimension, we take L\ as the invariant. 
This quantity is just the action or phase space area|l8||. 

dql2iT, where p and q are momentum and position 
coordinates, respectively. 

The work done by the particle as the piston extends 
from L to (1 + a)L, when the initial energy is i?o, is then 



= 1 - 



1 



(1 + a)^ 



(D3) 



The initial energy E^ is boltzmann distributed. Multi- 
plying the boltzmann factor by the density of states for a 
particle in one-dimension, dp/dE ^ 1/^/E, produces the 
following expression (apart from normalization factors): 



lim p{W,U) 

ts -^OQ 



1 exp f-/3M/ii±^ ) (D4) 
V a(2 + a] 



W 



which has the same functional dependence on W as (jDl|) 
when (3—1, apart from the overlap factor f{W). 

With the distribution (jPip . the Jarzynski relation is 
verified explicitly as follows. Using the integral 



s 



(D5) 



one obtains 

(e^^) = / P{W)e'^dW 
Jo 

^/a(2-ha)/7r 1 + a r(l/2) 



a 



^ + " [a(2+a)]i/: 



.(D6) 



1 + a = e-^-^ (D7) 
where AP = - In (1 -I- a) = - In (1 + AL/L). 
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FIG. 1: Definition of the forward (expansion) and reverse 
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FIG. 2: The structure of the overlap factor f{W) that modu- 
lates the exponential in the distribution function. This factor 
becomes a rapidly oscillating function in the limit of small 
piston velocities (i.e. a quasistatic process). 




FIG. 3: Probability distribution for the non-zero work values 
in the forward and reverse processes, when the piston moves 
fast. 



